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The confinement induced resonance provides an indispensable tool for the realization of the 
low- dimensional strongly interacting quantum system. Here, we investigate the confinement induced 
resonance in spin-orbit coupled cold atoms with Raman coupling. We find that the quasi-bound levels 
induced by the spin-orbit coupling and Raman coupling result in the Feshbach-type resonances. For 
sufficiently large Raman coupling, the bound states in one dimension exist only for sufficiently strong 
attractive interaction. Furthermore, the bound states in quasi-one dimension exist only for sufficient large 
ratio of the length scale of confinement to three dimensional s-wave scattering length. The Raman coupling 
substantially changes the confinement-induced resonance position. We give a proposal to realize 
confinement induced resonance through increasing Raman coupling strength in experiments. 

The experimental realizations of synthetic gauge field and the spin-orbit coupling (SOC) in neutral cold atoms 
provide a new arena to explore the exotic effects in cold atomic physics. For example, the SOC could bring 
about nontrivial ground states in Bose-Einstein Condensate (BEC), such as vortex or vortex lattice states, 
plane wave phase, standing wave phase 110 . The prominent effect induced by the SOC in fermions is that the SOC 
could enhance the low energy density of states, which results in the formation of two-body bound states and 
enhancement of the pairing gap 11,12 . In polarized fermion gas, the spin-orbit coupling modifies the finite tem- 
perature phase diagram 13 . The exotic dynamic effect, e.g. Zitterbewegung oscillation, appears in the spin-orbit 
coupled cold atomic gas 1415 . In the presence of SOC, the two-body scattering properties in three dimension have 
been investigated. It is shown that the SOC usually results in the mixed-partial-wave scattering 16 . For the low- 
energy scattering, the short range behaviors of wave function in three dimension can be modified by the SOC 1719 . 
The two-body scattering properties in qusi-two-dimensional confinement with pure Rashba spin-orbit coupling 
are investigated 20 . 

For the low-dimensional quantum gas, the two-body scattering properties can be affected greatly by the 
external confinement potential. For example, when the s-wave scattering length is comparable to the transverse 
confinement length (aja ± = 1/C with C = —((1/2) ~ 1.46), there exists a resonance, wherein the one-dimen- 
sional effective interaction constant diverges 21,22 . The similar scenario of the confinement induced resonance 
(CIR) also occurs in the quasi-two-dimensional case 23,24 . The confinement induced resonance has been observed 
through producing confinement induced molecules in quasi-one -dimensional Fermi gas 25 . It is also found 
experimentally that a single resonance splits into two resonances by introducing the anisotropic confinements 26 . 
The transversally anisotropic confinement alters the position of resonance by tuning the anisotropic ratio 27,28 . The 
confinement induced resonance which can be used to tune the interaction between atoms, provides a crucial 
ingredient to realize the strong interacting low-dimensional systems, such as Tonks-Girardeau gas 29,30 and 
possible Tomonaga-Luttinger liquid 31,32 . 

Some novel quantum states, for example, topological superfluidity, Majorana edge states or non-Abelian 
anyons could emerge in the low-dimension spin-orbit coupled quantum gas with Zeeman field 33-35 . In experi- 
ments, an effective Zeeman field in spin-orbit coupled atomic gas can be produced by two-photon Raman 
coupling 36 41 . The Raman coupling strength corresponds to the effective Zeeman field strength. The combination 
of the Raman coupling and SOC plays an essential role in the formation of the above novel quantum states. The 
effects induced by the Raman coupling and SOC are usually considered within BCS (Bardeen-Cooper-Schrieffer) 
mean-field framework 42,43 . It is known that the two-body interaction properties provide basis for understanding 
the many-body system. The studies on the effects of the Raman coupling on two-body problem may give some 
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insight into exotic quantum states. In addition, the confinement 
induced resonance provides the indispensable tool for the realization 
of the low-dimensional strongly interacting quantum gas. Further- 
more, how the Raman coupling and SOC affect confinement induced 
resonance is an inevitable question to clarify. In the present paper, we 
try to address the above questions by studying the two-body scatter- 
ing problem in one dimension and the confinement induced res- 
onance in the presence of the Raman coupling and SOC. 

Results 

The two-body scattering in the presence of spin-orbit coupling 
and Raman coupling. We consider the Hamiltonian of spin-orbit 
coupled cold atoms with Raman coupling 

H=H K +H 0 + V[x), 

K 2 

Hk— — -t — i«>- 



H 0 -- 



4m 

k 2 
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where the H K is the Hamiltonian in center of mass coordinate of two 
atoms, H 0 is the free Hamiltonian in relative coordinate. K and k are 
the total momentum and relative momentum of two atoms along x 
direction, respectively, f and f are the spin Pauli matrix and V(x) is 
the interaction between two particles, y denotes the SOC strength 
and Q is the two-photon Raman coupling strength between two 
Zeeman sublevels in experiment. The SOC strength is determined 
by y = 2nhsin(0/2)/X, where X is the Raman laser wave length, m is 
the mass of atom, 0 is the angle between two Raman beams. The 
above Hamiltonian is realized experimentally in fermion atomic gas 
of 40 K 38 . The above spin-orbit coupling is the mixture of equal 
Rashba and Dresselhaus type. Here, we choose the x axial direction 
as the direction of momentums. For simplicity, in the whole 
manuscript we denote K and k as the total and relative 
momentums, respectively, rather than K x and k x . The above 
Hamiltonian can be obtained from the Hamiltonian in Ref. 38 by 
applying a rotation in spin space (a x — » a z , a y — * a y and a z — * —(f). 
We take the natural units m = 1, H= 1 and y = 1 in this section. 

From Eq. ( 1), we know that the motion of center of mass is coupled 
to the relative motion through spins. In the following, we focus on 
the subspace of Hamiltonian with K = 0 (the effects of the non-zero 
total momentum are also discussed). In the spin basis 

[k) = (t4-|t)/\ / 2^i) = tt,|t2) = U,|t3) = (Ti + IT)/v / 5], the 

Raman coupling and SOC are transformed into 

M = 2y[fc/v / 2|5)((t 1 |-(t 2 |)+fe. C + Q/2>'(|t 1 )(t 1 |-|t 2 }(( 2 |)]. The 

interaction between cold atoms can be modeled by zero-range 
pseudo potential. Furthermore, we consider two identical spin- 1/2 
fermions. Hence, only the s-wave interaction in the singlet channel 
has contribution to two-body scattering. Therefore, the interaction 
matrix between two atoms takes the form as V(x) = \s)(s\ ®giD<>(x), 
where g lD is one-dimensional interaction constant. From the matrix 
V and M, we know that the spin channel 1 1 3 ) is decoupled from other 
channels and not affected by the interaction. Thus, the spin channel 
1 1 3 ) is dropped in the following, and the Hamiltonian H 0 is reduced to 
a 3 X 3 matrix. After diagnalizing H 0 , the eigenenerges are obtained 



as E 1 {k) = k 2 +2yJk 2 + (n/2yY, E 2 (k) = k\ E 3 {k) = k 2 



2yy k 2 + (fl/2y) (see the panel (a) of Fig. 1), respectively. When 

the total momentum K ¥^ 0, the spin channel | f 3 ) will couple with 
other states. Then, the H 0 is a 4 X 4 matrix and there are four 
continuous energy branches [see the panel (b) of Fig. 1]. 

The scattering problem can be solved through the Lipmann- 
Shwindger equation (for details see Metholds). The scattering with 
Raman coupling and SOC is intrinsically multi-channel scattering 



problem 44,45 . There exist different scattering thresholds for different 
energy branches. When the incident energy crosses the thresholds, 
some scattering channels are opened or closed. In certain scattering 
energy interval, there may exist several scattering channels scattering 
each other (see Fig. 1). The scattering amplitudes f m n (reflection 
amplitudes) make up a matrix of rank 1, where the subscript n (m) 
denotes the specific incident (reflecting) channel with a specific 
energy e. A single total amplitude / is obtained by diagnalizing the 
scattering matrix in every energy interval. 

It is known that in the usual case without Raman coupling and 
SOC, the reflection coefficient approaches one (total reflection) as 
the incident energy approaches the scattering threshold of e t h = 0- For 
a fixed incident energy e>0, the reflection also approaches total 
reflection as the interaction g 1D approaches infinity. For attractive 
interaction g 1D < 0, there always exists a bound state below the 
scattering threshold. 

In Figs. 2 and 3, we calculate the reflection coefficient ([/| 2 ) from 
the obtained scattering amplitude/. We can see that, compared with 
the usual case, the Raman coupling and SOC cause fundamental 
changes in the behaviors of the scattering amplitude at low energy. 
First of all, there exist scattering resonances in the parameter space 
because there exist quasi-bound states between the energy branches. 
The interaction matrix V can be rewritten in terms of the eigen-basis 
of H 0 . In addition to the interaction in the respective eigen-basis 
channel (diagonal part), there are also non-diagonal part coupling 
to different eigen-basis. The upper energy branches could support 
bound states near the incident energy. In addition, the bound states 
are coupled to the scattering states due to the non-diagonal inter- 
action term. Hence, it results in a Feshbach-type resonance scatter- 
ing 4648 . Secondly, when the collision energy approaches the lowest 
threshold, the reflection coefficient does not approach 1. So the 
incident wave could not be totally reflected in the presence of SOC 
and Raman coupling. Thirdly, the reflection may vanish under 
certain conditions. As shown in Fig. 2 and panel (a) of Fig. 3, 
as the incident energy approaches the threshold of e = 0 from 
below, the scattering amplitude becomes zero no matter how large 
the interaction is. This is because when the incident energy 
approaches the threshold of zero energy from below (e->0_), the 
middle energy branch E 2 contributes a infinitely large real number 
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denominator of the scattering amplitude (see Eq. (4) in Metholds). 
Next, as the interaction g 1D approaches infinity, contrary to the usual 
case, the reflection needs not mean a total reflection. Due to the 
coupling between singlet and triplet channels caused by the SOC, 
the incident wave can tunnel through other channels, even the inter- 
action strength is very strong in singlet channel. A similar scenario 
also appears in impurity scattering problems of multi-component 
coupled system 49 . 

Finally, the large Raman coupling changes the conditions of the 
existence of bound states. When the Raman coupling satisfies 0 < Q 
£ 2y, there always exists a bound state below the lowest threshold for 
attractive interaction g lD < 0 as the usual case. However, if Q. > 2y > 
0, there exists a bound state only when the attractive interaction is 
strong enough. This is because when the energy approaches the lowest 
£2), all the energy branches contribute a finite part 



threshold (e 
1 
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dk 



(fi/2) 2 



+ - 



e-E 2 +ii]k 2 +(n/2y e+^-£ifc 2 +(n/2) 



1+ 



1 



, to the denominator of the scattering 
e+h 1 -E 3 k 2 + (D./2) 2 1 

amplitude (see Eq. (4) in Metholds), rather than an infinitely large 
number as that in usual case without Raman coupling and SOC. The 
existence of the bound states is directly related to the density of states 
near the lowest scattering threshold. It is known that, in the usual case, 
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Figure 1 | The energy spectrum of the relative motion (the Raman parameter €1 = 4/Viri). (a): There are three energy branches (K = 0), the highest E\, 
the middle one E 2 and the lowest branches £ 3 . When the Raman coupling is weak (£2 < 2y), there exist two minimums on the lowest energy branch £ 3 (not 
shown in the Figure). When the Raman coupling is strong enough (£2 > 2y), the lowest threshold is — Q locating at the zero momentum of the 
lowest energy branch £ 3 . The A, B and C label three different scattering energy intervals [—4, 0], [0, 4] and [4, °°), respectively. The numbers of the 
scattering channels in different energy intervals are different, (b): Due to non- vanishing total momentum (K= 5y), the spin-state If J. + J. j) would couple 
with other spin states. So there would be four continuous energy branches. 



the density of states in one dimension near the scattering threshold 

(eg, = 0) behaves like ~ —= . There always exists a bound state bellow 

v e 

the scattering threshold as long as the interacting {g 1D < 0) is 
attractive (no matter how weak it is). When the Raman coupling is 
weak (Q £ 2y), the lowest energy branch has two minimums occur- 




e(7 2 /m) 



Figure 2 | The reflection coefficient f\ 2 as a function of the interaction 
gio and the scattering energy e when £2 = 4/7m (K = 0). /is the scattering 
(reflection) amplitude in one dimension. The corresponding three 
scattering energy intervals in panel (a) of Fig. 1 are also labeled here. 



ring at non-zero momentums. The density of states near the mini- 
mums also looks like cc —= (the e is measured with respect to the 

minimum of the continuous spectrum). 

Then the conditions of existence of bound states are the same as 
the usual case. However, if the Raman coupling is strong enough (fl 
> 2y), the lowest energy branch has only one minimum locating at 
the zero momentum (see the Fig. 1). Intuitively thinking, the density 

of states should be proportional to —=. However, at the zero 

momentum, we can see that the spin-orbit coupling does not play 
role at all in the energy spectrum (see Eq. (1)). The lowest point 
should correspond to the spin-triplet state, rather than the singlet 
state. Furthermore, we consider the interaction occurring only in the 
singlet channel. So the effective density of states near the lowest 
threshold which could contribute to the formation of bound state 
low. For example, there exists a suppression factor 

, in the denominator of the scattering amplitude 

k 2 + (n/2f) 5 

(see Eq. (4) in Metholds). The effective density of states near the 
lowest scattering threshold behaves like \fe approaching to zero as 
e— >0 (the e is measured with respect to the lowest scattering thresh- 
old), rather than like l/y/e blowing up. So, in this case, there exist 
bound states only when the interaction is strong enough. In the 
following, we will see that the modifications of effective density of 
states near the lowest threshold also have significant impacts on the 
existence of the bound states of quasi-one-dimensional system. 



very 
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Figure 3 | The reflection coefficients in the energy interval A of Fig. 1. (a): the refection coefficient corresponding to K = 0. (b): the refection 
coefficient corresponding to K = 5y. The Feshbach resonance peaks correspond to the quasi-bound states embedded between the continuous energy 
branches. The resonance position and width vary with the increase of the interaction strength. We can see that when the incident energy e approaches the 
threshold of the middle branch E 2 , the reflection vanishes. When the incident energy e approaches the lowest threshold the reflections are not total 
reflections. 



The modifications of scattering properties are robust even for non- 
zero total momentum 0 (see the panel (b) of Fig. 3). However, 
the position and width of Feshbach Resonance are usually modified 
by the non-zero total momentum. Fig. 4 shows the energy of bound 
states (with respect to the lowest scattering threshold) for attractive 
interaction (g 1D < 0) in the pure one-dimensional system. We can 
see that, with increase of total momentum K, the critical interaction 
magnitude \gm\ where the bound states begin to appear near the 
lowest threshold becomes larger and larger. From the above discus- 
sion, we know that for stronger Raman coupling fl, the required 
interaction strength \gio\ should be stronger in order to form a 
bound state. So the effects of the K 0 is similar to that of increasing 
Raman coupling. In the next section, we will see the conclusion that 
the increase of total momentum K amounts to increasing the Raman 
coupling Q is also valid for the quasi-one-dimensional system. 



psuedo-potential interaction between atoms 50,51 . g 3 £, = 4ph 2 a s /m 
and a) ± are the three dimensional interaction constant and 
frequency of confinement trap, respectively. a s and (i = mil are 
the s-wave scattering length and the reduced mass of two atoms, 
respectively. In this section, we take natural units as m = 1, H= 1 
and co± = 1. 

The one-dimensional effective interaction constant g lD can be 
obtained in terms of the three dimensional s-wave scattering length 
a s (see Methods), 



with 



gw- 



2h a, 



\ia\ l-(Ci + C 2 )a s /a_L ' 



(3) 



The confinement induced resonance. The one-dimensional effec- 
tive interaction constant is derived through investigating two-body 
problem of three dimension with confinement. The Raman coupling 
and SOC may change the condition of confinement induced 
resonance. After separating the motion of center of mass, the Hamil- 
tonian with confinement can be written as 



H T = H 0 +H L + V(r), 



Hj_ = 



h 1 



2^ + <>*) + ^ 



where H 0 is the free Hamiltonian of relative coordinate along x 
direction as above section, H ± is transverse confinement, 
V(r) = \s)(s\ ®g3D(>(~r)dr(r : ) is the three dimensional s-wave 



-1 
-1 



dt 



dt- 



ef'/ 2 1 

e et/2 



■1 



,-et/2 k 2 



dk- 



coshf yt^k 2 + (£l/2yf 



fc 2 + (n/2y) 2 

The above equation gives the connection between the one-dimen- 
sional effective interaction constant g 1D and the three dimensional s- 
wave scattering length. When the SOC vanishes (y = 0), C 2 becomes 
zero. Furthermore, if the scattering energy e— > 0, the constant C\ = C 
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Figure 4 | The energy of bound states for attractive interaction (g 1D < 0) in one dimension From the Fig. 4, we can see the energy of bound 
states is getting lower and lower with increase the interaction magnitude \gir>\ ■ The critical interaction magnitude is large for large total momentum K. 
So the effects of K are very similar to that of increasing Raman coupling. 



= —£(1/2) ~ 1.46, the resonance condition is reduced to Olshanii's 
result 21 . 

From panel (a) of Fig. 5, we can find that, with the increase of the 
Raman coupling, the resonance position a s /a ± is getting smaller. In 
addition, the resonance position inclines to be independent of SOC 
strength for sufficiently large Raman coupling (Q2>1). In fact, the 

resonance position a s j ai ~ 1 j \fl£l can be arbitrarily small by 

increasing the Raman parameter Si for fixed SOC parameter (y ~ 
1). It can be shown that, for a fixed spin-orbit coupling (y = 1), 
Q ~ VlQ diverges and C 2 is bounded as D, — » When the condition 
of confinement induced resonance is satisfied, from Eq.(3), the res- 
onance position a s j «i = 1 /(Ci + C2) ~ 1 / VlQ.. It means that, 

compared with the usual case, it is much easier to fulfill the condi- 
tions of confinement induced resonance in the presence of Raman 
coupling and SOC. For sufficiently large Raman coupling, the experi- 
mental observation of confinement induced resonance needs not 
resort to the usual magnetic Feshbach resonance techniques. 

For fixed SOC strength (y = 1), we show in Fig. 6 the bound state 
energies supported by the "closed" channels (the transversely excited 
modes) and the full Hamiltonian, respectively 22,27 . The confinement 
induced resonance can also be viewed as a Feshbach resonance as 
that in the usual case without Raman coupling and SOC. The res- 
onance condition is satisfied when the energy of the bound states in 
closed channel coincides with the scattering threshold of the ground 
transverse modes. In the meantime, the difference between the 
bound state energies in the closed channels and the full Hamil- 
tonian is 2hw±. 

However, based on the same reasons as that in the pure one- 
dimensional case, the large Raman coupling also changes the existing 
condition of the bound states in the quasi-one dimension. When the 
Raman coupling satisfies £1 £ 2y, there always exist the bound states 
irrespective of the sign of the s-wave scattering length a s , which is 
consistent with the usual case 25 (see the green and blue lines in the 
panel (a) of Fig. 6). However, if Q > 2y, the bound states exist only for 
sufficiently large a ± /a s (see red lines). As stated before, the reason is 



that the spin-orbit coupling and the Raman coupling change greatly 
the effective density of states near the lowest scattering threshold of 
energy spectrum. The spin-orbit coupling term yk (a\ — o\) does not 
play significant roles near zero momentum. So the states correspond- 
ing to the lowest and highest spectrum belong to spin-triplets, while 
the state in the middle energy branch is the spin-singlet. In addition, 
we consider the interactions occurring at the spin-singlet channel. So 
the bound states can be roughly viewed as the spin-singlet bound 
states which are mainly supported by the middle energy spectrum. 
Furthermore, the bound states exist below the lowest threshold (the 
threshold of triplet energy spectrum). So, the binding energy (relative 
to the middle energy branch) of bound states has to be larger than the 
differences between the thresholds of the middle and the lowest 
energy branches. In other words, the large Raman coupling Q pushes 
the energy of bound state low by amount of Q. at least. So, for strong 
Raman coupling, there exist bound states only when the interaction 
is strong enough. Then the bound states belong to the deeply binding 
states of the spin-singlet channel. In consequence, the position of 
confinement induced resonance is also modified accordingly. 

The panel (b) of Fig. 5 shows the dependence of the resonance 
position on the variation of the total momentum K. When the Raman 
coupling Q. = 0, the spin-orbit coupling terms can be removed by a 
gauge transformation. So the resonance position is not affected by the 
total momentum K (see the green line in the panel (b)). When the 
Raman coupling fl 7^ 0, the resonance position aja j_ becomes smal- 
ler and smaller with the increase of K. The panel (b) of Fig. 6 shows 
the energy of bound states for K 7^ 0 in the quasi-one-dimensional 
system. Similar to the effects of increasing the Raman coupling Q, the 
strong total momentum K also modifies the existing condition of 
bound states in quasi-one dimension. So the increase of total 
momentum K amounts to the increase of Raman coupling Q if the 
Raman coupling Q 0. 

Here we see what roles the confinement and the SOC (Raman 
coupling) play, respectively. The confinement freezes the transverse 
degree of atomic motion, and produces a quasi-one-dimensional 
system. Due to the enhancement of density of states, there is always 
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Q| 1 1 1 1 1 1 1 1 1 1 

01 23456789 

K(y/hmco±) 

Figure 5 | The resonance position aja ± as a function of the Raman parameter £1 and total momentum K (y = \Jhma)\}. (a): as shown in panel (a) of 
Fig. 5, when the Raman coupling is absent (fl = 0) and the incident energy approaches the lowest threshold (the minimum of £ 3 ), the resonance condition 
is also exactly the same as the case without SOC (aja ± = 1/(Q + C 2 ) = 1/C ~ 0.68). It is related to the fact that the constant gauge potential can be 
gauged away by applying a gauge transformation when the Raman coupling is absent. However, in the presence of non-zero Raman coupling, the 
resonance condition at the lowest threshold could never recover the usual case, (b) : the resonance position as a function of total momentum K. For non- 
vanishing Raman coupling (£2 # 0), the resonance position aja ± is getting smaller smaller with the increase of the total momentum K. 



two-body bound state in the quasi-one-dimensional system no mat- 
ter how weak the interaction is, as long as the interaction is attractive. 
After introducing the SOC and Raman coupling, the couplings could 
reduce, or even change fundamentally the effective density of states 
near the lowest threshold. The existence condition of bound states is 
also modified accordingly. The above two competing factors deter- 
mine the formations of bound states in the confined spin-orbit 
coupled system. 

Discussion 

The Raman coupling and SOC have been realized experimentally 
in 40 K atomic gas 38 . Two magnetic sublevels |f) = 1 9/2, 9/2) and \[) 
= 1 9/2, 7/2) are chosen as two spin 1/2 states. One can choose the 
experimental parameters y / m~2nh/ (ml) ~ 1. 28cm /s, io ± ~ 2n X 
17 kHz with a ± ~ 172 nm. The s-wave background scattering 
length a s ~ 170o 0 ~ 9 nm. Under the above conditions, the SOC 
strength y~ \/hmoj±, the ratio a s /a ± = 0.052. To satisfy the con- 
dition of confinement induced resonance, one can increase Raman 
coupling strength (Q ~ h X 3 MHz) by tuning the intensity of 
Raman beams. The binding energies of bound states near confine- 
ment induced resonance (Eg~2hw ± ) can be measured by using 
radio-frequency (rf) spectroscopy 52 . It is expected that there are 
two peaks in the radio-frequency photodissociation spectra. One 
peak locates at the atomic transition frequency (v 0 ) of an occupied 
state to another initially unoccupied state (e.g. 1 9/2, 7/2) — » 1 9/2, 5/ 
2)). The other peak at a non-zero detuning (S = — v 0 ~ 34 kHz) 



from the atomic transition corresponds to the dissociation of quasi- 
one-dimensional bound states 25 . 

In summary, we investigate the effects of SOC and the Raman 
coupling on the confinement induced resonance. The Raman coup- 
ling and spin-orbit coupling fundamentally change the interacting 
properties of atoms. We propose to realize the confinement induced 
resonance by increasing Raman coupling strength. Different from 
the usual way, such as utilizing Feshbach resonances to produce a 
large scattering length, our work gives a new way to realize the 
strongly interacting quasi-one-dimensional atomic gas with Raman 
coupling and spin-orbit coupling. Due to the exotic effects induced 
by Raman coupling and spin-orbit coupling, a lot of interesting 
many-body physical phenomena, e.g., the crossover of BCS-BEC 
superfluidity 53,54 , inhomogeneous Fulde-Ferrell-Larkin-Ovhinnikov 
(FFLO) state 55 , fermion pair breaking in the presence of external 
magnetic field 56 , need to be revised in the strong interacting quasi- 
one dimension atomic gas. 

Methods 

In this work, we investigate the two-body scattering problem in one-dimensional and 
the three dimensional cold atomic systems with confinement, respectively. Through 
comparing the scattering state in the one dimension with that in three dimension, one 
can get the effective one-dimensional constant (gio) m terms of the three dimensional 
s-wave scattering length (a s ). 

1 Two-body scattering in one dimension. In this section, we give the detailed 
calculation of the scattering amplitude. We consider the one-dimensional Hamiltonian 
of two atoms in the presence of the Raman coupling and the spin-orbit coupling 
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cq -2 



-4 - 



(b) 



K = 0 
- - if = 5\/fimuj± 




■ 

a±/a 



1 2 

s 



Figure 6 | The bound state energies in the closed channels and the full Hamiltonian, respectively (with SOC strength y= \/hmco±). (a): The 
bound state energies are measured with respect to the lowest thresholds. The thin lines denote the bound state energies supported by the closed channels, 
and the thick lines are the bound state energies in the full Hamiltonian. With the increase of Raman coupling £2, the resonance position aja ± (the 
reciprocal of a ± la s ) is getting smaller and smaller, which is consistent with that in panel (a) of Fig. 5. (b): the energy of bound states for various total 
momentums K. We can see that with the increase of the total momentum K, the modifications of existing conditions of bound states are very similar to 
that of increasing the Raman coupling. Note here the green line in panel (b) corresponds to the blue line in panel (a). 



H — H K + H k , 



The interaction between atoms can be modeled by zero-range pseudo potential. 
Therefore the interaction matrix between two atoms takes the form as 



4m 2m • 



-«?). 



V(x) 



fg lD &(x) 0 0 0\ 

0 0 0 0 

0 0 0 0 

\ 0 o o o/ 



h*=-+^-(os-«^)+?(»;+»o+v. 

m m • 2 ' 

When the total momentum K — 0, the Hamiltonian of center of mass (H K ) is zero. In 
the spin basis [(T|-1T) / V^.(TT)»(U).(tl + it) / v^], the Schrodinger equation 
in the relative coordinate is 

[H 0 + V(x)]V = EV, 



k 2 

H 0 = — +M, 



with 



0 
0 



where g 1D is the one -dimensional interaction constant. From the interaction matrix V 

and matrix M, one can get that the spin channel (t! + 4T) j is decoupled from 

other channels in the case of K — 0. So we will drop it in the following. The 
Hamiltonian is reduced a 3 X 3 matrix. After diagnalizing the H 0 , the eigenstates (in 
coordinate space) and eigenenerges are obtained 



•Mao 



Vk 2 + & 2 - 



\ 



A J 



*3Jfc(*) = 



M = 2y 



0 


k 

Vi 


k 

v5 


k 

•Ji 


A 


0 


k 


0 


-A 


0 


0 


0 




We introduce A — — in the above equation. 



s/2k 



\ Vk 2 + A 2 + A J 
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corresponding to eigensnergies Ei (k) — k 2 + 2y\J k 2 + A 2 , E 2 (k) — A; 2 , 

E^(k) — k 2 — 2y\J k 2 + A 2 , respectively (in the natural units of m — 1,^ — 1 andy — 1). 
The energy spectrum is shown in panel (a) of Fig. 1. 

The scattering problem can be solved through the Lipmann-Shwindger equation 

>F(x) =%(*)+ \dx'G(ejc,x')V(x') x l'(x'), 



where G(e,x,x') = ( x 



x' ) is the free Green's function. In general, 



e-E m + it]\ 

the Green's function is a 3 X 3 matrix. In the case of zero range interaction the 
Lipmann-Shwindger equation takes the following form 



/G„(*)\ 



G 21 (x) 



1-£idG„(u) 

where is the first component of incident state T 0 (x) and 



1 f°° 



dkxe'" 



_J k 2 

e-E 1 (k) + h]2(k 2 + A 2 ) 



1 



e-E 2 (k) + ink 2 + A 2 

1 k 2 

' e-E,(k) + it]2(k 2 + A 2 ) 

27tJ_„ [e-^ + i^^ + A 2 

e-/c 2 /c 2 
" (e-/t 2 + iij) 2 -4y 2 (/t 2 + A 2 )*: 2 + A 2 



1 f°° 

G2i(x) = — J dfcxe* 



e-EiW + ifj 4(fc 2 + A 2 ) 

1 - v^/tA 1 

e— E 2 (k) + iij 2(k 2 + A 2 ) + €—E 3 (k) + iij 



1 



-V2Ak 



4(k 2 + A 2 ) 

i r°° 

: — \ dkxe ikx , 
2jt J_ QCl [e-fc 2 + if/2(fc 2 + A 2 ) 

V2[Afc(e-/c 2 )+2yA:(J: 2 + A 2 )] 1 
(e-k 2 + i n f-4y 2 (k 2 + A 2 ) 2(fc 2 + A 2 ) 



(4) 



G 21 (x)-- 



\f2Ae 



V2 [A (e - cc 2 ) + 2y (fz 2 + A 2 )] e*> 1,1 
4(* 2 - a 2 )( a 2 + A 2 ) 



4(a 2 + A 2 ) 

^[A(e-z 2 )+2y( a 2 i + A 2 )]e i °M 
4(a 2_ a 2 )(o; 2 + A 2 ) 



G 31 (x) = 



s/2Ae 



^2[A{ 



-2y(x 2 + A 2 )]e""M 



4(«| + / 
V2[A(6- 



4(a 2 -« 2 )(«f+A 2 ) 



-2y(a 2 + A 2 )]e''"'W 



4( a 2 -a 2 )(o : 2 + A 2 ) 



with X, = y e + 2y 2 - 2y \Je + y 2 + A 2 , x 2 = v'i and a 3 = \j e + 2y 2 + 2y \ < 
For x > 0, the above equations is expressed in another form 

(Gu(x)\ 

\ G 3lW/ 



where/i = — ^J 2 



iA 



->fs- 



For a given incident energy, 



2y, 



there is possibility that there exist several states scattering each other (see Fig. 1). From 
the Eq. (4), we can explicitly obtain the scattering amplitude for a given incident state. 
The scattering amplitude makes up to a scattering matrix of rank 1, it can be reduced 
to one single total amplitude/by diagnalizing the scattering matrix. According to the 
number of the scattering states, there are several distinct cases. 

The incident energy is above the highest scattering threshold Q<€ [within energy 
interval C (see Fig. 1) ]. When incident energy satisfies Q. < e, there will be three eigen- 
states scattering each other forwardly (all a t , a 2 an d a 3 are real number), which locate 
at the right semi-axis of relative momentum (k > 0) (see Fig. 1). The resulting total 
scattering amplitude is 

_ gw [*u m + *u m + *u m] 

f ~ l-gwG n (0) ' 



The incident energy is above the zero-energy, but below the highest threshold Q 
(0 <e<Q. within energy interval B). When incident energy lies in 0 <e<f2, the 
scattering state on the highest branch is closed (oil becomes imaginary number). The 
remainder two eigen-states scattered each other forwardly. The scattering amplitude 



I-SidGh(O) 



dk x e' k 



e-E^ + nj 
V2k(yk 2 +A 2 -A) , 



4(fc 2 + A 2 ) 

-V2kA 1 
2(k 2 + A 2 ) + <E-E i (k) + in 

\/2*:(V*: 2 + A 2 + A) 
4(/c 2 + A 2 ) 



e-E 2 (k) + in 



■ dkxe ikx 
2n J 



-\/2A/c 



e-k 2 + ii]2(k 2 + A 2 ) 
V2[Ak{e-k 2 )-2yk(k 2 + A 2 )] 1 



(e-A: 2 + ^) 2 -4y 2 (yc 2 +A 2 ) 2(k 2 +A 2 ) 



with n — » 0+. The above integrals can be obtained by using the residue theorem 

AV*'* (e-orJ)aie*"H (e-f) 2 )^'"^ 



G n (x) = i 



2=< 2 (a 2 + A 2 ) 2{a\-c l 2 ){a 2 + A 2 ) 2(v 2 -a 2 ) (a 2 + A 2 ) 



The scattering energy satisfies — Q<e<0 (within energy interval A). When the scat- 
tering energy £ lies in the regime — Q < e < 0, there is only one scattering state on the 
lowest energy branch (both a x and y. 2 become imaginary number). The scattering 
amplitude is 



/ = 



l-gi D G n (0) 



2 The confinement induced resonance. The scattering problem in three dimension 
with confinement can be solved as follow. We assume the confinement is strong 
enough that only the transverse harmonic ground state is occupied. The incident 
energy with respect to the lowest threshold and the ground state energy of transverse 
harmonic oscillator should be lower than the transverse exited state energy. The 
lowest excited state which can be coupled to ground states by the s-wave interaction is 
<j) i(y)<j> i(z) with energy 2ho)± 22 - 57 . Then, the incident energy of scattering states should 
€<2hcQj_ — 2 (in natural units of m — l,h— landco x = 1). Using the Lipmann- 
Shwindger equation, three dimension scattering wave function can be obtained 

^(r), D = Uy)U^«^)+ \dr'G w (e,rj')V(r')V w {r'), 



where c^„(f) = (-^n2"n\a ± ) 



-!/2„ 



,2 / 2 <H„(t/. 



"±),a ± = 



--\[2\ 
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function and the length scale of transverse Harmonic oscillator. T 0 (x) is one 
dimensional incident wave function along x direction as above section. The 



interaction Vis the pseudo -potential V(r) 



471/7 a 

I s) (s I ® S(r)d r (r). Substituting 



the interaction V, the wave function in three dimension becomes 

/(GdM)iA 

fm(.x,y,z)=<l>o(y)<l>o(z)'¥Q(x)+g3DF\ (G 3 D(r)) 21 , 

V(G3DM) 3i y 

where F — lim r _>od r [rf^P^fr)] > ^^(r) is first component of three dimensional wave 
function. Focusing on y — z — 0 and the resulting qusi-one dimensional wave 
function is 



¥(*) 



: ^0(0)^0(0) 



l-fiD(Qu(0) u ) f 



/(G 3D M,0)) u \ 

(G 3D (*,0,0)) 21 
V(G 3D (x,0,0)) 3 J 



(5) 



where (030(0)!!),. = /fm r ^. 0 5 r [r(G3 D (r) n )] is the regular part of the Green's function 
matrix element Gn at origin. It is known that when r —> 0, the three dimensional 
Green's function in confined system also diverges as that in homogeneous space 
(G 3 £>(f*)ii)r-»-o °c — lMTrr 58 . In order to obtain the regular part, we need to subtract the 
singular part (— l/4nr) from G 3 o(^)n near the origin. For the qusi-one dimensional 
scattering problem, we only need to know the long-ranged asymptotic behavior of 
three dimensional Green's function and the regular part of G 3 d(^»^ — 0, z = 0) n near 
x — 0. The Green's function can be decomposed as two terms: 



G 3D (e 3 x,)' = 0,z = 0) 

- E 

— G 3 ni - 



x,0,0 



I 



0,0,0 



e-H 0 + H ± + ni 
e-[ni + « 2 + l-E m (fc)]+i»; 



where 



and 



I )t 2 + A 2 ' 

It can be shown that the /2M is not singular as x — » 0, while the J£t(ac) d° es diverge 
Xifjc) (x — l/4nx. Substracting the diverging part, 



Zlr(0) = - 



Jik r- 

8^J„ 



dl 



So the regular part of Green's function (G 3D (0)„) r = G u (0)|^(0)| 4 + / lr (0) + x 2 (0). 
For jc>1, the Green's function G 3d (jc>1,0,0)~G 3d1 (x,0,0) = (0) | 4 G(x), G(x) is 
the Green's function in pure one domension. Comparing Eq. (4) with (5), we can get 
the one-dimensional effective interaction constant in term of the three dimensional s- 
wave interaction constant 



gw = 



l-«lD(Zlr(0)+X2(0))' 



AH 2 na s 



Introducing ^ lr (0) — \flC\ ^ Sn, x 2 (0) — \/2C 2 j 87t and using g^ D 
a± — \/2, m — 2/.1 and A — — , the above equation is reduced to Eq. (3) in the main 



giD- 



2h 2 a t 



tia\ l-(C 1 + C 2 )a s /a± ' 



with 



- 1 r 



dt 



' <f>l 2 j_ 
.(e'-l)Vt t 3 / 2 . 

0l 2 



V2 



1 r 00 ef'/ 2 r°° 
Jo dt ~l\- 

yt\/t 2 + (a/2y) : 



rffcr* 2 '/ 2 



J /c 2 + (fl/2y) 2 



G 3D i = 2^ 



e-[i-E m (k)]+iti 



and 



G3D2 — 



(■GO 

_ V j te -[».+"! + i-i!»,(t)-e]< 
fc(0),C W»,(°)<W 0 )<* I ^aX^ I 0)1 



here x P m j fe ) are the eigenstates of H 0 which is same as that in the above section. The 
summation fc ^ excludes the term [«j — 0, « 2 ~ 0, fc, m], the identity 

1 



n 



r 

Jo 



4«i,n 2 fc, 

dte~ ni (n > 0) has been used in the above Equation. The first part G 3D1 



is 



related to the scattering channel and long ranged, while the second part G302 is short 
ranged corresponding to influences of virtual transition of other closed channels. 
Redefining the reference point of energy e— >€— 1, it is easy to see that the first term 
G 3D1 is just the one dimensional Green's function G{x) up to a factor |0o(O)| 4 = l/2n. 
In the following, we will see the short-ranged part of matrix element {G^rxiixjiu is 
related to confinement induced resonance (CIR). 



ing the summation over the index | 
(G 3 D2(*))ii takes the form 



, the second part of the matrix element 



(G3m(x)) n -- 



1 r» <?l* r 



8?t 2 
A 



k 2 + A- 
Z1 + Z2. 



+ cosh\ytVk 2 + A 2 



k 1 



k 2 + A 2 



where 



ZiW = 



1 r°° i^ 1 / 2 c°° 
~\ dt—\ dke- k '" 2 + ,k ', 
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